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SUMMARY 

A numerical technique has been developed at Langley Research Center which 
generates linear perturbation models from nonlinear aircraft vehicle simulations. 
The technique is very general and can be applied to simulations of any system 
that is described by nonlinear differential equations. The computer program 
used to generate these models is discussed with emphasis placed on generation 
of the Jacobian matrices, calculation of the coefficients needed for solving 
the perturbation model, and generation of the solution of the linear differen- 
tial equations. Included in the paper is an example application of the tech- 
nique to a nonlinear model of the NASA Terminal Configured Vehicle. 


INTRODUCTION 


Linearized models of physical systems, when used either in conjunction with 
nonlinear simulations of the physical systems or in an independent mode, offer 
the research engineer many insights to his problem. Obviously, a nonlinear sim- 
ulation will be a more valid model over a much larger range of the system's 
state variables, but with the present state of the art of mathematics and sys- 
tems design techniques, linear models offer many advantages. 

A linear model used to represent the system over some limited region has a 
known analytical solution which can be programmed on a digital computer and does 
not require standard numerical integration techniques (ref. 1). This property 
can result in a savings in both computation time requirements and computer mem- 
ory allocations for the simulation of the system. Many available computer algo- 
rithms have been written which will identify the eigenvalues and eigenvectors 
of a linear model (ref. 2). This gives the reser'cher a quick look at the char- 
acteristic modes of the system. Once these modes have been identified, steps 
can be taken to eliminate undesirable characteristics by adding a feedback con- 
trol system. At present, many books and articles have been written on the sub- 
ject of linear feedback controls, and many computer programs are available which 
will provide such features as root placement and the solution to both the time- 
varying and steady-state optimal regulator problems (refs. 3 to 7) . 


This report will describe a computer program which was designed to obtain 
linear models about a nominal state and control vector from nonlinear real-time 
aircraft simulations. The program is very general in design and may be applied 
to any system that is described by a set of nonlinear differential equations 
about any trajectory in state space. The program uses various Lagrange interpo- 
lation formulas to obtain both the state and control Jacobian matrices. Once 
they are obtained, the linear differential equations are integrated by using the 
local linearization technique described in reference 1. Eigenvectors and eigen- 
values are calculated using standard computer routines that are available at 
Langley Research Center. 
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n x n dimensional state Jacobian matrix, A(t) - — 


X-Xq 


element of A(t) located at intersection of ith row and jth column 


n x k dimensional control Jacobian matrix, B(t) 


*4 !■ 


mean aerodynamic chord, meters 


F vector of total aerodynamic forces 

f(x,u) n dimensional vector of general, nonlinear time-varying functions of 
state vector x and control vector u 

f^( ) ith component of f 

g acceleration due to gravity, meters per second 

h integration interval step size, seconds 


i 

' i 


h(t) n dimensional vector whose elements are residual higher order terms 

from Taylor series expansion 

I n x n identity matrix 

I XX» I YY' I ZZ» I XZ moments of inertia, kilograms-meters^ 

J u n x k Jacobian matrix of f with respect to u 

J x n x n Jacobian matrix of f with respect to x 


k constant equal to number of elements in control vector u 

L/D lift-drag ratio 

8.(5 ( ) coefficients used in Lagrange interpolation formulas for approximation 
of fi 

£(J { ) coefficients used in Lagrange interpolation formulas for approximation 
of Sfj/Sxj 

M vector of total aerodynamic moments 

M^v maximum operating Mach rmmhor 
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n constant equal to number ot elements tn state vector x 

P n v n dimensional matrix which i a Pad * 1 approximation to 

A* 1 ( *Ah . I)R 

**••0 par lod of characteristic mode, seconds 

Pta roll tat*, radians per second 

Q n ' k matrix which is a Pad** apptoximat ion to A“‘(e* h - Ah - l)» 

% pitch tat* , radians per second 

r^ yaw rat*, radians per second 

S aim lari ty t tans for mat ton matrix 

T thrust, newtons (MtC. in computer -genet ated tables) 

t time, an independent vat table, seconds 

tf final time, seconds 

tj start ins) t ime, second* 

t)/j time to damp to one-half amplitude, seconds 

1 2 time to double amplitude, second* 

% longitudinal translation velocity, meteis pet second 

u(t) k dimensional vectot whose elements are contiol variables of system 

V c cal ibt ated a it speed, knots 

v max maximum opeiatlng airspeed, knots 

lateral translation velocity, motets pet second 
w^ vet t leal translation veUx Ity, motets pet sevx'ud 

x ( t ) n dimensional vectot whose elements ate state vat tables of system 

xJ nominal state vectot with jth element possibly dittetent ttixn its 

nominal value 

•S small pet t urtvtt ion of vatiable away ttv>m its luxninal value (OKI, iu 

cvwput et -genet ated t ables) 

Vj, atieton position, deg tees 

elevatoi position, degt ees 
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6j constant which is amount jth element of nominal state or control vec- 

tor was perturbed 

rudder position, degrees 

6 8 stabilator position, degrees 

<5sp,L flight spoiler position, left side, degrees 

^sp,R flight spoiler position, right side, degrees 

9 pitch attitude, radians (THETA in computer-generated tables) 

5 DH damping coefficient for Dutch roll mode of aircraft 

£p damping coefficient for phugoid mode of aircraft 

£sp damping coefficient for short period mode of aircraft 

t variable of integration 

Trs time constant of roll subsidence mode 

tsd time constant of spiral divergence mode 

<J> roll attitude, radians (PHI in computer -generated tables) 

yaw attitude, radians (PSI in computer-generated tables) 

Subscripts: 

o nominal values of variables 

R rotor 

A dot over a variable indicates a time (DOT in computer-generated tables). 

PROBLEM DESCRIPTION 

Aircraft simulated on Langley's real-time simulation system are described 
by a set of nonlinear simultaneous differential equations of the form 


X ( t ) * f [x ( t ) ,u(t) ,t] 


( 1 ) 


where t represents time, x(t) is an n dimensional time-varying state vec- 
tor, u(t) is a k dimensional time-varying control vector, and f is an n 
dimensional vector of general nonlinear functions. As shown by reference 3, if 
UQ(t) is a given input (control) to the system described by equation (1) and 
Xo(t) is a known solution of the state differential equation, one can find 

approximations to neighboring solutions for small deviations from the initial 
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state and input vectors by using a linear state differential equation. Assume 
that XqU) satisfies 


*o(t) - f[xo(t) , Uq ( t) ,t] 


(tj i t i t f ) 


and that the system is operated close to nominal conditions. Therefore, one 
can write 


u(t) - u c (t) + 6u < t ) 


x(t) • XqU) + 6x(t) 




J 


(ti i t l t f ) 


( 2 ) 


Substituting equations (2) into the state differential equation (eq. (1)) and 
expanding in a Taylor series about (Xq < t ) , U q ( t ) ) yields 


x Q (t) + <5x(t) - f^Xjjtt) ,Uo(t) ,t] + J x [x 0 <t) ,u D (t) ,t 1 <$ x ( t ) 

+ J u [x 0 (t),u 0 (t),t] 5G(t) + h (t) 

where J x and J u are the Jacobian matrices of f with resp ct to x 
respectively. They are given by 

9f I 


(3) 

and ii , 


A £ 


3x 


B 5 


x-*o 

u-u rt 


?f 

3u 


x»x c 

u«u. 


The term h(t) is the sum of higher order terms from the Taylor expansion and 
should be "small" with respect to «5x and <Su. Neglecting h(t), 6x and iSu 
approximately satisfy the "linear" equation 


<5x(t) - A (t) flx (t) + B (t) <5u ( t ) 


(4) 


iSx ■ A6x + Bi5u 


(5) 




which is called the linearized state equation. For the particular applications 
of interest, only the time invariant system is considered in which A(t) and 
B(t) are constant matrices A and B. Therefore, the linearized state equation 
can be written as 
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NUMERICAL LINEARIZATION TECHNIQUE 
Computation of the A and B Matrices 

Now by using equation (1) and by assuming that each component (fj (x,u,t) 
for i “ l,n) is continuously differentiable m times and can be evaluated m 
times, the partials required for the A and B matrices can be approximated 
by using the Lagrange interpolation formulas (refs. 8 and 9). The components 
of f(x,u,t) are approximated by 
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(i - 1 #n) 


Due to notation complexity, this formula is explained by an example that uses 
the three-point Lagrange formula. First, 


" ( x ot» x 02 » * * *» x j» • • *» x o n ) 


which is the nominal state vector with the jth element allowed to vary from 
its nominal value while all the elements remain fixed. For the three-point 
formula (m ■ 3) , these vectors are 


X 1 ^ " ^ x 0 ’j» x 02 * * • • > x o j ~ ^j* • • •» x o n j 

x 2^ ” ^ x 0]* x 02' * * • ,x Oj » • • • ,x Op^ 


x-jj ■ 


^ x oi' x 02* • • ** x Oj + ^j* • • ‘' x o n ) 


(Xj - Xo^Txj - (x 0j + 6j)' 

^ ( X Oj - 6 j) - X Oj |"( x Oj - 6 j) - ( X Oj + 6 j)] 


(Xj - XOj) ( x j - "Oj - h ) 


^2 (xj) ■ 


Xj - (x 0j - «j)] [xj - (x Qj + « j) I 
’xoj - ( x 0 j - « j)i r*o - (x 0j + «j) 


(Xj - x 0j + «j)(xj - Xo. - 6j) 


1 I I I 
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[xj - (Xoj - fij)] (Xj - Xoj) 
l3( ’' 3 ’ ’ [(»Oj ♦ «j) - («0) - «j)] [(«Oj * S l) - «0)‘ 

(Xj - Xoj ♦ S j)(Xj - «Oj) 


In order to compute equation (6) must be differentiated with rospect 

3x-; - - 

: x-xq 
U«Uq 

to Xj and the resulting equation evaluated at (x 0 ,Uq) as follows: 


a i j S 


^ x 3 x=x Q k«l 


yi ^lc(Xj) fifXfcjjUQ) 


Again using the three-point formula as an example yields 




2(xj - Xq^ - fij 


(Xj) 


- 2 ( x j “ x Oj) 


^3 (Xj) 


2 (Xj - Xq^ + fij 


which, when evaluated at the three values of Xj and summed according to equa- 
tion (7) , results in 

3 f i } r , ■ i 

~ * ^|_- f i( x 1 j f«o) + f i(*3 j r5o)J ,8) 

X 3 x-Xq 3 
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An equivalent result for the five-point differentiation formula is 

3fi 1 j. _ 

7~ « rTT f i( S l j » 5 o) * 8f i(x 2 j « i5 o) + 8fi(x43,u 0 ) - f i (x 5 3,u 0 ) 

3x4 _ _ 126s L J 


and that for the seven-point differentiation formula is 


kjlS-So 60 6 j 


— [-f i (x 1 3,u 0 ) + 9f i (x 2 3,u 0 ) - 45fi(x33,u 0 ) + 45fi(x 5 3,u 0 ) 


- 9fi(x 6 3,u 0 ) + fi(x 7 3,u 0 ) 


The computation of the B matrix is identical to that of the A matrix 
except that x is held constant and u is varied. 


Use of the Perturbation Model 

Once the A and B matrices have been determined, the perturbation model 
defined by equation (5) is ready for use by the researcher. The two most com- 
mon uses of this model are to use it in place of a nonlinear simulation in stud- 
ies that will be very limited in their area of operation and to determine the 
eigenvalues of the model about the defined state trajectory. These eigenvalues 
in aerodynamic problems identify the basic modes of the aircraft; these are 
Dutch roll, short period, phugoid, spiral divergence, and roll subsidence. In 
this program, the eigenvalues of the perturbation model (characteristic roots 
of the A matrix) are determined by a standard Langley library routine. 

If the researcher desires to use the perturbation model in place of, or 
to compare with, his nonlinear model, it will be necessary to integrate equa- 
tion (5). The solution to equation (5) is 


6x(t) * e At 6x(o) 


•f 


e A(t-T) Bu(T) dT 


where the solution to the nonlinear system would be approximated by 
x(t) * x(o) + 6x(t) 

As shown by reference 1, a discrete approximation to equation (11) with local 
truncation error good to O(h^) is given by 

6x| {+ i * e^ixfc + Pfiufc + Q6uj< (1 
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where 


P - A* 1 (e^ - I) B n 3) 


Q • A" 2 (e Ah - Ah - I)B (14) 


5u - 5 u(5_) 

<Su k 5 ( 15 ) 

h 

and as before 

x k+l ■ x k + '■ Sx k+1 (16) 

Equations (12) to (16) are solved by the program’s integration subroutine. 

PROGRA USAGE AND LIMITATIONS 

For normal application the following should be followed by the user 
(fig. 1): 

(1) Trim the nonlinear aircraft model about the desired trajectory to 
obtain the nominal state vector x Q and the nominal control vector Up. The 
trim algorithm used at Langley Research Center for most real-time simulations 
is described in reference 10. 


XNCM 


(2) Compute the A matrix by using subroutine JACMAT (appendix A) with 


'o 


and M * N a n (number of states) 


(3) Reset the states to their trim values. 


(4) Compute the B matrix by using JACMAT with XNOM = U 0 , M * k (number 
of controls), and N * n. 

(5) Reset controls to their trim values. 

(6) If eigenvalues are require’, the user must call a subroutine which gen- 
erates eigenvalues. For the applications presented, subroutine REQR, a part of 
the Langley computer mathema: * ca 1 library, was used. 

(7) If integration of the linear system is required, call subroutine COEFF 
(appendix A) with NDIMA - n and NCOLB » k for calculating the coefficients 

of $Xk, <5 Uk, and v '' u k* 

(8) Obtain response of the linear system to a predetermined input sequence 
uk by calling subroutine INTEGRT (appendix A). 
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When applying this technique to general nonlinear simulations, certain 
potential problem areas should be mentioned. First, all implicit loops in the 
nonlinear equations must be broken by substituting variables and by reformu- 
lating the equations. If this is not possible, an iterative technique may pos- 
sibly be used to determine the approximate perturbed steady-state forces and 
moments. Second, the magnitudes of the perturbations used for the state and 
control variables need be chosen with care because if the perturbations are too 
small or too large, the derived linear model will not be a good approximation 
to the nonlinear analysis. The method used to choose perturbation magnitudes 
for the NASA Terminal Configured Vehicle (TCV) example is described in appen- 
dix B. And third, even though the linearization technique can be applied about 
any nominal state trajectory, th? results are more meaningful when the vehicle 
is trimmed and the nominal trajectory is stable. 


PROGRAM APPLICATION 


As an example of the results obtained from a standard application, the TCV 
airplane, a Boeing 737-100, was chosen. The desired outputs of this applica- 
tion were (1) linearized models of the B-737 about various trim conditions, 

(2) identification of the basic modes of the aircraft (eigenvalues) at these 
trim conditions, and (3) time-history comparisons of the linear and nonlinear 
models for predetermined inputs. 

The desired linearized models were of the form 
6x * A6x + B6u 


where the state vector was chosen to be in the body-axis system. The el 
of the body-axis state vector xp are 


u b 

w b 

% 

0 


*b 


vb 


Pb 


c b 

$ 

* 


10 



and the elements of the control vector are 


Tables I to V are example outputs and show the nominal state and control vec- 
tors, the body-axis A and B matrices, the eigenvalues, and the correspond- 
ing eigenvectors. 

Linear models defined in other axis systems can be derived from this model 
by means of a similarity transformation S where 

x b - Sx D 


x b * Sx D 

with x D being the desired state vector defined in the new axis system and S 
being time invariant. Substituting into our linear model 


yields 


6x L Afxfc + B^u 


S6xp * ASfxp + Bfu 


6x n * S -1 AS6xn + S _1 B6u 


as our linear model in the desired axis system. 


To further show the usefulness of these linear models as a simulation veri- 
fication and validation tool for the various flight conditions shown in table VI, 
a comparison of independent Boeing data (unpublished) and the linear models gen- 
erated from tho nonlinear simulation is shown in tables VII ano VIII. A review 
of these tables will show that good agreement exists between the simulation 
models and the independent data in almost all cases, excluding the spiral diver- 
gence mode. However, major disagreeme ts do exist in the short period mode of 
condition V with the aft center of gravity and in the phugoid mode of condi- 
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tion VII with the forward center of gravity. The researcher should now take 
steps to resolve the reasons for the differences in these cases by first reveri- 
fying the implementation of the nonlinear simulations' aerodynamics data in 
these areas and by trying to obtain other independent data such as flight data. 

Additional insights into the system being simulated can be gained by com- 
paring the linear models generated by each Lagrange interpolation formula. A 
general indication of the linearity of the simulation about the nominal trajec- 
tory is obtained, as well as an indication of sensitive modes and parameters 
(nonlinearities) of the simulation. For example, a comparison of the models 
obtained for the maximum speed case (table VI, condition V) showed that the 
lateral and the short period modes were approximately linear, but for the phu- 
goid mode, the damping ratio varied by 36 percent and the natural frequency by 
3 percent. This information implies that the linear models obtained would not 
be suitable for studies requiring precise knowledge of the phugoid mode. 

The numerical linearization technique has also been successfully applied 
to nonlinear simulations of other aircraft. Linear models of a fighter air- 
craft, a general aviation aircraft, a standard rotorcraft, and the rotor systems 
research aircraft (RSRA) developed by NASA and the U.S. Army (ref. 11) have been 
obtained about various trim conditions. The standard procedure as described was 
used in all cases except that of the RSRA aircraft which required procedural 
modifications since the nonlinear simulation included a dynamic rotor model 
which was continuously rotating during the linearization. Figure 2 outlines 
the iterative technique used to obtain the linear models for this vehicle. 
Basically, the approach taken was to allow integration of the rotor dynamics. 
However, a steady-state condition had to be obtained after each perturbation 
of a state or control before numerically calculating the Jacobians. Averaging 
the forces and moments over a number of rotor revolutions and at various points 
during each revolution is also done to enhance the credibility of the linear 
model obtained. 


CONCLUDING REMARKS 

The numerical linearization technique described in this paper has been 
successfully applied to nonlinear simulations of various aircraft. At this 
writing, linear models of the NASA Terminal Configured Vehicle, a fighter air- 
craft, a general aviation aircraft, a standard rotorcraft, and the RSRA have 
been obtained about various trim conditions. Linear models of aircraft with 
stability augmentation systems have also been obtained by augmenting the state 
vector with the associated automatic control system states and by proceeding 
in the manner described in the paper. 

A modification of the technique for application to simulations of rotor- 
craft with dynamic rotor models has also been developed and described. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 18, 1978 
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APPENDIX A 


DESCRIPTION AND LISTINGS OF SUBROUTINES JACMAT, CORFF , AND INTEGRT 

The major portion of the linear analysis package consists of three subrou- 
tines, JACMAT, COEFF, and INTEGRT. 


Subroutine JACMAT 

Purpose : Calculates the Jacobian of the nonlinear system about a nominal point 

in vector space. This corresponds to the A matrix when the input argument 
is the state vector and the B matrix when the input argument is the control 
vector . 

Use : CALL JACMAT (XNOM,F,FVAL , JACOBN , DELTA, N ,M,EOM, 1 FNTS ,MAXROW,MAXOOL) where 

XNOM An N dimensional input vector; this vector contains the nominal 

values of the independent variables about which the Jacobian is 
calculated 

F An N dimensional output vector; duiing computation of the partial 

derivatives, it contains the values of the dependent variables 


FVAL 


JACOBN 


An N " M ' IPNTS dimensional a: lay used in calculating the partial 
derivatives; FVAt.(l,J,K) is the Ith component of F evaluated 
at trie Kth change in the Jth component of XNOM 

An N > M dimensional output array which is the Jacobian matrix eval- 
uated at XNOM; that is. 


JACOBN ( t , J) 


Of i 
Jx i 


x-XNOM 


DELTA 

N 

M 

ROM 


An N dimensional input vectoi to step sires; PF.LTA(l) is the 
increment for XN0M(1) 

An integer input specifying the number of equations 

An integer input specifying the number ot independent vaiiables 

A user-supplied subroutine which calculates the values of F used 
in computing FVAL; ROM is a subroutine in the parameter list of 
JACMAT r the statement 

RXTRRNAL ROM 

must be included in the calling p root. am of JACMAT; the calling 
statement for EOM is 


1 t 









APl'KNDIX A 

CALL BOM (N ( M,XNOM,P) 

vt.<-r» N, M, end XNOM ere input* end F is the output 

IPNTS An integer input which specifies the Intel point ion formula to tx* 

used t 

IN ITS - 2, three-point formula 
1NFTS • 4, five-point formula 
TNPT? « ti, seven-point formula 

MAX NOW An integer input specifying the maximum number of equations to be 

used 

MAXCOL An integer input specifying the maximum numbei of independent vaii- 
ablee to tx> used 

The listing of subroutine J ACM AT is as follows. 
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APPENDIX A 


Subroutine COEFF 

Purpose : Computes the coefficients (e Ah , P, and Q) required for calculation 

of the discrete approximation to the solution of perturbation model. 

Use : CALL COEFF (A,NDIMA,B,NCOLB,H ,EAH r P,Q,W,MAXDIMA,MAXCOLB) , where 

A An NDIMA * NDIMA dimensional input array? this array is the A 

matrix of the perturbation model 

NDIMA An integer input specifying the dimension of A 

B An NDIMA * NCOLB dimensional input array; this array is the B 

matrix of the perturbation model 

NCOLB An integer input specifying the number of columns of B 

H Length of the integration interval 

EAH An NDIMA * NDIMA dimensional output array which approximates e Ah 

P An NDIMA * NCOLB dimensional output array which approximates 

A -1 (e Ah - i) B 

Q An NDIMA * NCOLB dimensional output array which approximates 

A -2 (c Ah - A h - I ) B 


W An NDIMA x NDIMA dimensional working space array 

MAX DIMA An integer input specifying the maximum dimension of A 

MAXOOLB An integer input specifying the maximum number of columns of B 

The listing of subroutine COEFF is as follows: 
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APPENDIX A 
Subroutine INTEGRT 

Purpose : Generates solutions to the linear differential equations obtained from 

the nonlinear simulation. 

Use : CALL INTEGRT (N,X,X0,X0D0TH,L,U,U0,UD0T,EAH,P,Q,W1 ,MAXN,MAXL) , where 

N An integer input specifying the number of states being used; 

N * MAXN 

X A MAXN-dimensional input/output vector which contains the values of 

the states in its first N locations; on input, X contains the 
past values of the states, and on output, it contains the current 
values of the states 

XO A MAXN-dimensional input vector that contains the values of the 

states at which the A and B matriceswere calculated in its 
first N locations; that is, XO contains X D 

XODOTH A MAXN-dimensiona? input vector which contains C*H in its first N 
locations where C is the value of f(x 0 ,u 0 ) and H is the same 
as in COEFF 

L An integer input specifying the number of controls being used; 

L S MAXL 

U A MAXL-dimensional input vector which contains the current values 

of the controls in its first L locations 

UO A MAXL-dimensional input vector that contains the values of the con- 

trols at which the A and B matrices were calculated in its 
first L locations; that is, UO contains U Q 

UDOT A MAXL-dimensional input vector which contains the time derivatives 

of the controls in its first L locations 

EAH A MAXN x MAXN-dimensional input array; this is the same EAH as in 

COEFF 

P A MAXN x MAXL-dimensional input array; this is the same P as in 

COEFF 

Q A MAXN x MAXL-dimensional input array; this is the same Q as in 

COEFF 

W1 A MAXN-dimensional vector used for working space 

MAXN An integer input specifying the maximum number of states 

MAXL An integer input specifying the maximum number of controls 


The listing of subroutine INTEGRT is as follows: 


21 


mf mm '+*#****<? 


| aw. 


4*1 "J4 | ■ » w'J * 


V I 






•1 ' j j 


1 i 


I 


\ 



V 



APPENDIX A 


•h (V <*) <« IT <0 


cpeoH(\in4ini0i-it9OH(vin4ip«Mi(roHiMn«ir 

^4 I ^^« r H M ^ l _t^, a ^ r ^ ( xiiV(V(V<viVi(V(ViV(Vi*>(*)(*)r)r r )r*i 


►—*—*—*— 

1- 

k- 

k- 

kkk 

k- k- k- 

k- k- 

k- 

k- k- 

h H 

k- k- 

k- 

h- ►- 

k- k- 

k- 

h h K 

k- k- k- 

arc arc 

X 

I 

X 

uccmri 

ct a. tr 

x x x x 

X T 

X 

a a 

X X 

X 

a x x x x x x 

CI900 

© 

© 

© 

m5(5l50(5t5l50(J(50l5(5ai5ffi00(S 

0(5 15(5(50(20 

UJ UJ UJ UJ 

UJ 

UI 

UJ 

LJliiUJUiUiliiU/U^UiUJUiU. UJ 

* J UJ UJ UJ UJ UJ UI 

UI UJ UJ UJ UI bJ U! U.' 

K k h K 

k- 

k- 

k- 

k— k— k— 

k- k- k- 

k— k- 

k- 

k- 1- 

k- k- 

k- k- 

k- 

k- k- 

k- | — 

►— 

kkkk 

k- k- k- 

zzzz 

z 

z 

z 

z z z z z z 

Z Z 

z z z z z z z 

z z z 

Z z 

z 

z z z z z z z 

M M M M 

M 

44 

44 

W H H 

M M H 

44 44 

H 


44 44 

M H 

H 

44 44 

44 44 

N 

44 44 44 44 

44 44 44 

_ «J 

c 

♦ 















X 41 

X 44 

k- 
















z ►* 

Z 

x 















• z 

o 

44 

N 














Z UJ 

►4 

X 















x x 

k- 

— J 

k- 

• 













X UJ 

3 

UJ 

o 














z u. 

_l 

o 

a 

-1 













• u. 

O 

* 

X 

X 

• 












*4 44 Z 

V) 

I 


< 

44H. 












X o o 


X 


X 

z 












• 1-4 

'w 

UJ 

o 


X 












3 UJ ►» 

X 


k- 

3 

< 

• 











• IX 

k- 

u 


» 

X 

X 



• 








a k- 3 



z 


*- 




3 








* 9 


. — 

o 

-J 

o 

k- 











X O UJ 

4 

» — 1 

44 

X 

X 

—1 



k- 








X k- 

o 

♦ 

k- 

X 


3 






• 





UJ _J 

UJ 


3 

z 


V) 



3 



X 


z 



• Z X 

> 



4 • 

z 

UJ 



V, 





X 

• 


k- O 44 

_J 

X 

o 

o Z 

X 

IK 



UI 



k- 


< 

44 


C k— 

o 

-J 

V! 

X X 

4 




a 



_J 


X 

3 


C k Z 

IS) 

UJ 


X 

z 

UJ 






3 





3 3 U! 


c 

bJ 

♦ z 


X 



u 



VI 


11 

k- 


• ja 

IT 


X 


X 

k- 



X 



u 



_ 


C) O UJ 

44 


k- 

44 O 

• 




k- 



X 



3 


3 V) U 


> 


#4 4 

44*4 

_l 









IT 


• U. 


OC 

• 

♦ 44 

z 

_J 



_l 



UJ 


44 

bj 


3 UJ 

— 


z 

SC JX 

X 



_l 



X 


> 

X 


• I c 

e 

z 

bJ 

— X 

X 

C- 



X 



»- 


X 



3 k- 

3 

UJ 

X 

X < 

z 




u 





X 

UJ 


• UJ 

* 

> 

k- 

-I z 

44 

Q 






-J 


X 

I 


I «)X 

o 

44 


UJ * 

44 44 

z 



a 



-1 


X 

k- 


k- UJ k- 

X 

(9 

4 

o z 

3 Z 

X 



z 



X 

z 

X 



o *- • 

*■4 


X 

X 

• X 




X 



u 



_l 


O X k- 

U 

VI 

• 

N X 

44 X 








N 

N 

-J 


o z in 


44 

h- 

X 

-1 z 

c 



c 



o 



X 


x —a. 

♦ 


o 

44 44 

x — 

X 



3 



z 

44 


44 U 


• X H 


z 

o 

44 a 

x x 



►X 


H 


X 

ro 


4^ 

o c u. 

3 

c 

X 

♦ 4 

X k- 

1 


< X4> 

1 





<4 

• c 

44 

x a 

_1 

44 


SC 4 

- o 



o 


o 



V 

> 

X Z 

3 

• a • 

UJ 

k- 

♦ 

~ z 

o o 

X 


X 

3 

3 


X 

X 

X 

• X 

• 

x a — 

c 

X 


X X 

3 C 







-J 

X 

X 

X 

_J 3 

• X 3 

• 

3 

#4 

X 

• X 

H 


1 

H 

1 


UJ 

X 

X 

X 

X • 

Z • 


o 

if 

z 

44 • 







O o 

X 

X 

UJ 3 

x a 

~ U> x 


UJ 

■w 

>- 4 

— J 4* 

X 


»44v 

3 

*«. 


• a 

— sc 

X 

• _! 

J z * 

k Z 

♦ 


k— 

CD Z 

x 

_J 


M 

_! 

»-X 


X 



> UJ 

>- 

a hu. 


-1 

o 

X 

x — 

UJ 


"»■* 

UJ 



X H 

N H 

n 

X o 

mux 

© k- 

X 

X 

o 

Z X 

z >- 

c 


X 

c 

3 


UJ 



X * 

X 

UJ 3 N 

-J 

44 

3 

UJ Z 

— X 







4- 

44 44 

44 

x a 

— — ct 

»- O 

bu 

k- 

* 

> 4- 

k- a 



N 


N 


44 

X IV 

IT 

X 

r x x 

2 #k 

a 

z 

c 

►4 X 

o a 

Uj 



UI 



Uj — 

44 4- 


X bJ 

~ — X 

►4 IT O 

• 

UJ 


© x 

C X 

k- 



k- 



k- > 

>» > 

> 

w k- 

> > *- 

3 C 

X 

a 

♦ 

UJ 3 SC 

X 



X 



X X 

X X 

X 

V) X 

XXV) 

uj in x 


U' 


V) 


-J 

z 

^"4 

_J 

-1 


-j a x x 

X 

X -J 

x x a 

Z 

N 

u 

-4 

k4 Z 


3 

• 


3 

• -- 


3 X 

x a 

a 

O 3 

X x o 

4- V) Z 


u 

x 

o 


O 



a 



<J X 

X X 

X 

k- Cl 

X X k- 

k- H O 

k- 

44 


44 


-J 

N 

X 


N ~ 


_) sc 

X X 

X 

X _l 

XXX 

3 I ►" 

O 

c 

3 

— V) 


X 

H 


X 

►-» 3 


X 



X X 

X 

O k- k- 

o 


-J 

3 Z 


<_> 



o 



u 



u 


a x 

X 

VI 

UJ 

• UJ 



*- 



a 





-1 

_J 

OC 3 

_l 

44 

c 

X Z 











_l 

_J 

3 O 

UJ 

X 

* 

4» 44 



o 



o 





X 

X 

V) U. 

c 

k- 

a 

u c 



c 



o 





u 

<J 


' (V 


— (V 

22 


mi W****! ai.;,J,-...v .J -j.-i J. ., < 



! 


<#WfSes*asfc®» 




' • l 

r 1 

t . 

i i 



r 


.4 


APPENDIX A 


O t- ® O' © •»• <Nj 

m n m n 4 4 4 


PI4lT<0K8 9o 

4 4 4 4 4 4 4 * fl 


ia xama 

(9 (9 (5 (9 O O (9 

llIU W Ul IiMiJ 111 


sasssxKX 

(5019819000 
111 II 1 111 111 li : 111 III 111 


ZZZ?ZZZ zzzzzzzz 


ORIGINAL PAGE IS 
OF POOR QUALITY 



.. h 


’ 'l 


' s 





f f 


.4 



\ 


APPENDIX B 


I METHOD OF CHOOSING PERTURBATION MAGNITUDES 

I 

j The magnitudes of the perturbations used for the state variables in the 

TCV example were chosen as a function of aircraft states in which the aeronauti- 

i i cal engineer would have some intuitive feel as to their desired range of varia- 

; \ tion. For this example, the variables used were total velocity V^, angle of 

attack a, angle of sideslip 8, roll attitude 0, pitch attitude 9, and yaw 

ii i attitude 0. The values of the variations in the body-axis state variables 

| are given by 

I u b 

i Auj;, = AViji wfc Aa - Vfc cos a A8 

V T 


w b 

Awb = Av t — + ujj Aa - vjj sin a A8 
V T 


i 


i 


Aq b = Aa sin 0 cos 9 + rb A0 + Pb sin 0 A9 


v b 

Avb * Av t — + V T cos 6 A8 

V T 


t 


Apb = -Aa sin 9 - a cos 9 A8 



Arb = Aa cos 0 cos 9 - qb A0 + pb cos 0 A9 

where 

„ g tan 0 
a - 

V T 


24 


Aa = 


g A0 

V T cos^ 0 


— AVm 
V T 






\ 





! 





AB - 0.1/57.3 rad 


AO • 1./57.3 rad 
At « 1./57.3 rad 


At - 1./57.3 rad 

It should be noted that all variables except t are at these trim values. The 
magnitude of t was not allowed to be less than 2/57.3 rad sr that a nonzero 
value for Aq b would be calculated. 

The variations in the control variables were 1 percent of the tc‘al range 
of each control variable. 
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TABLE I.- VALUES OF HIE STATES, CONTROLS, AND STATE DERIVATIVES 
AT THE POINT AT WHICH THE PERTURBATION MODEL 
WAS GENERATED 

["Aircraft weight, 36 287.4 kg; altitude, 457.2 m; airspeed, 

I 63.09 m/sec; flap deflection, 40°; landing gear down 

STATE STATE DERIVATIVE CONTROL 


UB 

3 

63.04486 

UBDOT 

3 

0.00001 

ENG 

3 

35339.07894 

WB 

m 

2.8632 

WBDOT 

= 

-0.00019 

STAB 

— 

8.56032 

QB 

= 

0. 

QBDOT 

3 

-0.00008 

DELR 

= 

0. 

THETA 

- 

-0.00698 

THETADOT 

B 

0. 

DELE 

= 

2.66833 

VB 

B 

0. 

VBDOT 

3 

0. 

DELA 

= 

0. 

PB 

= 

0. 

PBDOT 

3 

0. 

SPL 

= 

0. 

R3 

XX 

0. 

RBDOT 

3 

0. 

SPR 

= 

0. 

PHI 

= 

0. 

PH I DOT 

K 

0. 




PSI 

S 

0. 

PS I DOT 

3 

0. 








! ■ 

TABLE II.- THE A HATRl v j 

l ' 

t 

[Aircraft wight, 36 387.4 Kg> altituda, 467.2 «» airapaad, j 

j 63.09 at/saei flap daflaction, 40°j landing gaar down j » 


UBDOT 

at/aac/aac 

-0.03781 

0.11139 

-2.86186 -9.80664 

0. 

0. 

0. 

0. 

0. 

\! * 4 

f . 

WBDOT 

at/aac/aac 

-0.28511 -0.72225 

63.01435 

0.07896 

0. 

0. 

0. 

0. 

0. 

1 

i 

li \ 

< X 

QBDOT 

rad/aac/aac 

-0.00055 -0.01971 

-0.50162 -0.00032 

0. 

0. 

0. 

0. 

0. s 

t 

i 


THETA DOT 

rad/aac 

0. 

0. 

1 .0 

0. 

0. 

0. 

0. 

0. 

*• !' ; 


VBDOT 

at/aac/aac 

0. 

0. 

0. 

0. 

PBD0T 

rad/aac/aac 

0. 

0. 

0. 

0. 

It BOOT 

rad/aac/aac 

0. 

0. 

0. 

0. 

PHIDOT 

rad/aac 

0. 

0. 

0. 

0. 

PSID0T 

rad/aac 

0. 

0. 

0. 

0. 


UB 

at/aac 

NB 

at/aac 

08 

rad/aac 

THETA 

rad 


• ! 
* 



0.14763 

3.39181 

-62.73394 

9.80633 

0. 

0.06835 

-1.88160 

0.99852 

0.00003 

0. 

0.01064 

-0.14972 

-0.14574 

-0.00452 

0. 

0. 

1.0 

-0.0069755 

0. 

0. 

0. 

0. 

1.0 

0. 

0. 

VB 

PB 

RB 

PHI 

PS1 

m/aac 

rad/aac 

rad/aac 

rad 

rad 





TABLE III.- THE B MATRIX 


l Aircraft weight, 36 287.4 kg; altitude, 457.2 m; airspeed, 

f | 63.09 m/sec; flap deflection, 40°; landing gear down 

» ; 

i UBDOT 0.00003 0.00458 0. 0.00220 0. -0.00252 -0.00252 

] av'sec/sec 

*! I 

’ j WBDOT 0. -0.10089 0. -0.04851 0. 0.02570 0.02570 

i m/sec/sec 

h QBDOT 0.00001 -0.04098 0. -0.01972 0. 0.00097 0.00097 

j ; r ad/sec/sec 


THETADOT 

rad/sec 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

: 

VBDOT 

m/sec/sec 

0. 

0. 

0.04154 

0. 

0.00026 

-0.00383 

0.00383 

! 

PBDOT 

rad/sec/sec 

0. 

0. 

0.01084 

0. 

0.01946 

0.01156 -0.01 1 56 

i 

‘ j , 

1 

RBDOT 

rad/sec/sec 

0. 

0. 

-0.01102 

0. 

0.00163 

0.00141 -0.00141 

PHIDOT 

rad/sec 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

i 

: \ 

PSIDOT 

rad/sec 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

i i 


THRUST 

N 

STAB 

deg 

DELR 

deg 

DELE 

deg 

DELA 

deg 

SPL 

deg 

SPR 

deg 

: 

1 



1 29 



* 4 - 


TV ITU mi 


\ J 


i 2 


■r .• -I 


'/ H \% ' 


i i 




TABLE IV.- EIGENVALUES OF SYSTEM 

Aircraft weight, 36 287.4 kg; altitude, 457.2 m; airspeed, 
63.09 m/sec; flap deflection, 40°; landing gear down 


EIGENVALUES 


fc. 


TIME DAMPING UNDAMPED NATURAL PERIOD t}/ 2 
CONSTANT RATIO FREQUENCY 


-.2016E+01 + 0. *1 . 4960E+00 

-. 5940E-02 + 0. *1 . 1 684E+03 

0 . + 0 . *1 
-.1635E-01 + . 1778E+00*! 

-.1635E-01 + -. 1 778E+00*I 
-.7636E-01 + .1138E+01*I 

7636E-01 + -.11 38E+01 *1 
-.6145E+00 + . 1 110E+01 *1 

-.614SE+00 + -.1110E+01*I 
ERROR RETURN FROM REQR = 0 


30 


.9161E-01 .1 785E+00 


.6694E-01 .1 141E+01 


•4845E+00 . 1 268E+01 


.3438E+00 

.1167E+03 

. 3. 35E+02 . 4238E+02 

. 5520E+01 . 9077E+01 

. 5663E+01 .11 28E+01 


1 >1 


. - a 



* 

\ 


\'\ y 


> i i 


•* [• -i ‘ » ■’ \m 


■ i .? 

4 - i 


TABLE V.- THE A MATRIX EIGENVECTORS 

weight, 36 287,4 kgi altitude, 457.2 m; airspeed, 63.09 a/secj^ 
flap deflection, 40°; landing gear down 


. 92388E-1 3 

-.7721 7E-11 

0 . 

-. 33790E+00 

-. 9341 9E+00 

-. 66972E-1 2 

-. 27383E-1 2 

-. 1 6303E-01 

. 1 5728E-01 

-. 31828E-11 

. 30459E-1 1 

0. 

. 2221 5E-01 

.112166+00 

-. 81 036E-1 1 

. 50657E-1 1 

-. 87648E+00 

. 48086E+00 

. 25839E-13 

-. 1 0793E-1 5 

0. 

-. 48066E-03 

-. 8901 6E-C3 

-.54591 E-l 3 

-. 29485E-1 3 

-. 30577E-02 

-. 44365E-02 

-.128156-13 

.1 81 71 E-l 3 

0. 

-. 471 89E-02 

. 31 381 E-02 

-. 22585E-1 3 

. 49478E-1 1 

-.1 891 9E-02 

. 38035E-02 

-.99308E+00 

. 288966+00 

0. 

. 31 229E-1 5 

-. 26995E-1 4 

-.655586+00 

-.754996+00 

-.237616-15 

-. 1 2345E-1 4 

-. 1 0499E+00 

-.18419E-03 

0. 

-.10451 E-l 6 

. 24668E-1 6 

. 891 03E-02 

.468586-02 

-.108776-16 

.591106-16 

-.65566E-02 

. 5681 SE-02 

0. 

-. 33435E-1 6 

-. 70731 E-l 6 

-. 291 04E-02 

.287636-02 

-.1 07356-16 

-. 753116-18 

• 52050E-0I 

. 37683E-01 

0 . 

. 76340E-1 6 

.8571 4E-16 

. 35568E-02 

-.808496-02 

.455086-16 

-. 1 4068E-1 6 

.3251 9E-02 

-.95658E+00 

• 10000E+01 

-.37741 E-l 5 

. 22281 E-l 5 

. 26866E-02 

. 23768E-02 

. 3581 0E-1 7 

. 7691 7E-1 7 









TABLE VI.- BASIC FLIGHT CONDITIONS 


Condition 

Description 

Weight, 

*9 

Center of gravity 



Flap 


Altitude, 

Mach 

V C' 

knots 

Forward 

Aft 








I 

Approach 

*40 016.3 

9.10c 

0.31c 

40° 

Down 


0 

0.185 

122 

II 

Melding 15 







Up 

Up 

1 

254 

.382 

230 

III 

Maximum dynamic 
pressure 











3 

962 

.831 

440 

IV 

Cl imb c 











3 

048 

.622 

340 

V 

Maximum speed 
and Mach 
number 











6 

096 

.900 

403 

VI 

Cruise 











6 

096 

.735 

330 

VII 

v max/**max 

corner 











6 

102 

.840 

350 

VIII 

Maximum altitude 
cruise 











10 

058 

.742 

250 

| 

IX 

Lightweight 

Vmax/^max 

corner 

d 31 746.0 



0.33c 





7 

102 

.789 

330 


a For weight of 40 823 kg 

I xx - 508 432 kg-m^ 

Iyy * 1 186 340 kg-m 2 

I zz » 1 626 981 kg-m 2 

I xz “ 105 754 kg-m 2 

^Chosen as being 10 pe-cent above maximum L/D speed. 
c Minimum cost climb. 
d For weight of 31 751 kg. 




TABLE VII.- CHARACTERISTIC MODES FROM BOEING DATA 



Center of 

Short period 

Phugoid 

Dutch roll 












v^onQ x t ion 













^SP 

p sec 

t l/2 

Cp 

p sec 

*1/2 

£dr 

p sec 

*1/2 

i 

0.1 

0.418 

5.39 

1.29 

0.084 

34 

45 

0.057 

4.81 

9.31 


0.3 

0.583 

8.58 

1.32 

0.130 

48 

40 

0.047 

5.07 

11.99 

ii 

0.1 

0,361 

2.86 

0.82 

0.030 

62 

225 

0*117 

. 

3.34 

3.12 


0.3 

0.490 

4.22 



0. '3 

0.028 

71 

273 

0.109 

3.58 

3.59 

hi 

0.1 

0.426 

1 .96 

0.46 

0.640 

118 

16 

0.125 

1.85 

1 .62 


0.3 

0.743 

4.67 

0.46 

0.263 

44 

18 

4 

0.121 

2.00 

1 .82 

IV 

0.1 

0.354 

1.97 

0.57 

0.065 

109 

185 

0.102 

2.49 

2.68 


0.3 

0.499 

3.04 

0.58 

0.086 

147 

188 

0.093 

2.70 

3.20 

V 

0.1 

0.419 

2.40 

0.57 

tl/2 « 5. 


0.134 

2.12 

1 .73 






t 2 - 

24. 9 a 




0.3 

0.924 

12.70 

0.58 

t 1/2 - 4. 

4? 

0.132 

2.30 

_ 

1.91 






to * 

17. l a 



] 

VI 

0.1 

— 

0.333 

1 .86 

0.58 

0.134 

' 

155 

' 

127 

- * 

0.089 2.32 

. f 

2.86 


0.3 

0.484 

2.89 

0.59 

ti/2 - 80.3? 

0.086 2.51 

3.19 






t 2 ■ 358. 

i a 




VII 

0.1 

0.366 

2.28 

0.64 

■ 

0.925 

L_ 

273 

12 

0.097 2.21 

_ . i. ..... 

2.49 


0.3 

0.645 

5.00 

0.65 

0.482 

60 

12 0.092 

t 

U f - -» 

2.37 

2.83 

VIII 

0.1 

0.285 

2.45 

1 .01 

0.103 

86 

92 0.083 

1 

3.08 

4.10 


0.3 

0.371 

3.76 

1 .04 

0.099 

87 

96 0.079 

3.30 

4.60 

IX 

0.1 

0.401 

2.18 

0.55 

tl/2 

- 17. 

1 

3? 

0.109 | 2.10 

2.11 






t 2 > 

7.2 a 





0.3 

0.739 

5.54 

0.56 

0.451 

50 

11 10. 097 

2.25 

2.56 

_ _ 











a Complex conjugate pair splits into two simple poles. 
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0.3 0.547 3.46 0.584 t-[ / 2 * 6 . 35 ; 

tl> 2 " 7.9® 

a Ccmplex conjugate pair splits into two simple poles 
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TABLE VIII.- Concluded 


Center of 

Condition gravity 


Spiral divergence Roll subsidence 

Root T SD| t 1/2 or t 2 Root t RS *1/2 or t 2 
3 312.5 

-1 429.0 
-37.0 

5 

6 

-111 .11 


5; 

5; 


0.003 


56. 

,77 

37. 

,59 

n . 

,94 

73. 

,53 

27. 

,32 


223.6 


77.0 


36.5 
36.5 
22.9 
96.7 
25 . 4t 
26. 0< 
49.8’ 


0.521 

0.361 

0.515 

0.357 

0.456 

0.316 

0.455 

0.315 

0.341 

0.236 


2.621 0.382 


2.09 0.479 


2.542 0.3931 


-2.538 0.394 


2.21 8 0.451 


2.217 0.451 


-1 .710 


2.638 0.379 


0.266 


.263 










































































